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ABSTRACT 



Aims. We investigate whether stellar dust sources i.e. asymptotic giant branch (AGB) stars and supernovae (SNe) can 
account for dust detected in 5 < 2 < 6.5 quasars (QSOs). 

Methods. We calculate the required dust yields per AGB star and per SN using the dust masses of QSOs inferred from 
their millimeter emission and stellar masses approximated as the difference between the dynamical and the H2 gas 
masses of these objects. 

Results. We find that AGB stars are not efficient enough to form dust in the majority of the z > 5 QSOs, whereas SNe 
may be able to account for dust in some QSOs. However, they require very high dust yields even for a top-heavy initial 
mass function. 

Conclusions. This suggests additional non-stellar dust formation mechanism e.g. significant dust grain growth in the 
interstellar medium of at least three out of nine z > 5 QSOs. SNe (but not AGB stars) may deliver enough heavy 
elements to fuel this growth. 
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1. Introduction 

Studies of the extragalactic background light have re- 
vealed that roughly half of the energy emitted in the 
Universe apart from t he CMB is reprocessed by dust (e.g. 
iHauser fc Dwekl |2001[ ). Thus, understanding the physical 
processes responsible for the formation of dust throughout 
cosmic time has important cosmological implications. 

Dust can either be formed by asym ptotic giant branc h 

(AGB) stars (even at low metallicities; ISloan et"al]l2009t) . 

or supernovae (SNe). Alternatively, the bulk of the dust 
. mass accumulation may occur in the interstellar medium 
' (ISM) on dust seeds produced by stars. This process 
' can successfully explain gas depletions in the Milky 

Way (iDraine fc SalDetei][l979l: iDwek fc ScalolfToSOt iDraind 
I1990I I2009D. alo ng with the dust 

(iMatsuura et all l2009l) and a z - 
IDwek et al.lT2007l) . 

Theoretical works have shown that an AGB 
star and a SN p r oduce up to ~ 4 x 10~ ^ Mq 
(|Morgan fc Edmunds! 120031: iFerrarotti fc Gaill I2006D a nd 



The dust in the Milky Way was predominantly formed 
by evolved stars with only a minor SN contribution 
(iGehrzl I1989D . but individual SNe may form significant 
amounts of dust. Submill imeter observation s of th e SN 
remnant s Cassiopeia A (iDunne et al. _ 20031 2009() and 
Kepler ([Morgan et all l2003t iGomez et all l2009t) have re- 



vealed as much as ~ 1 Mq of freshly formed dust, but 



Krause et al 



these results are controversial ([Pwekl [200^ 
200 4 | Gomez et al.ll2005l:lAYilson fc Batrlall2005HBlair et al. 
20071: ISibthorpe et al.l2009HBarlow et al. 120101) . Dust yields 



10-3-10-2 MfT 



(iGreen et al. 2004 Borkowski et al. 2006 


: Sugerman et al. 


120061: lErcolano et al. 120071 iMeikle et al. 


20071 Rho et al. 


!2008'.'2009: Kotak ct al.'i2009': 'Lcc et al.| 2009jSakon et al. 


12009: Sandstrom ct al. 2009: Wesson et all 20091). 



masses 
6.42 



of the LMC 
quasar (QSO; 



1.32 M^ lflbdini fc Ferrarall2001t iNozawa et al.ll2003[) of 

dust, respectively. However, for the case of SN dust, only 
< O.lM fTi of the dust actuall y surv iv es in the associated 
shocks (jBianchi fc Schneider! l2007t ICherchneff fc Dwekl 
I2OIOD . 



The situation is even more complex at high redshifts. 
iDwek et"aI1 (|2007t ) claimed t hat only SNe ca i i prod uce dust 
on timescales < 1 Gyr, but IValiante et al.l (j2009t ) showed 
that AGB stars dominate dust production over SNe as 
early as 150-500 Myr aft er the onset of star formation. 
iMichatowski et all (j2010bl ) concluded that in three out of 
six 4 < z < 5 submillimeter galaxies, only SNe are efficient 
enough to form dust provided that they have high dust 
yields. This would then be suggestive of a significant dust 
growth in the ISM and/or a top-heavy initial mass function 
(IMF). 
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Table 1. Dust, gas and dynamical masses of 2: > 5 QSOs 













Mdyn sin^ i 




Tdust 


No. 


QSO 


z 


(10'Mq) 


(IO^OMo) 


(10"A/o) 


(K) 


1 


J0338+0021 


5.03 


7.1+0.6 


2.2 " 


3.0 " 


31 


45.6 " 


2 


J0840+5624 


5.85 


4.7+0.9 


2.5 


24.2 ^ 


53 




3 


J0927+2001 


5.77 


7.2+1.1 


1.8 


11.8 


25 


51.1 


4 


J1044-0125 


5.74 


2.7+0.6 


0.7 


0.8 


26 




5 


J 1048+4637 


6.23 


4.3+0.6 


1.0 


4.5 


23 


<40 " 


6 


J1148+5251 


6.42 


5.9+0.7 


1.6 " 


4.5 " 


27 


55.0 ^ 


7 


J1335+3533 


5.93 


3.4+0.7 


1.8 ^ 


3.1 


53 




8 


J 1425+3254 


5.85 


3.3+0.7 


2.0 * 


15.6 


60 




9 


J2054-0005 


6.06 


3.4+0.8 


1.2 


4.2 


35 





Notes. The s ample inclu des QSOs detected in their dust continuum and CO l ine emission ("Carilli et al.' '2000", '200l!. l2007l: 
Bertoldi et al... 2003a. b: P etric et aljr20 03: Priddey et al. 2003, 2008; Wal ter et all [20 03. 2004; Robson et al. 2004; Beele n et all 
20061 : iMaioliiio et al.ll2007l : iRiechers et aLll2007i : IWan g ct al. 2007, 2008a, b. 20101 : 1 Wu et al.ll2009f l. We calculated dust masses from 
the detection of 1200 pm emission using Eq. 5 of [Micha lowski ct al, QOOy) assuming /3 = 1.3. The errors reflect the statistical 
uncertainties only. For three QSOs we assumed observationally inferred dust temperatures (the last column). For the rest we 
ad opted the average o f the se estimates Tdust = 50 K. iMaiolino et al.l (|2007f ). '''' IWang et al] (|2010l '). ''^^ [Walter et all l[2004f ). 
f'*' IWang et all l[2008bD . iRobson et aP ([20041 ). (•'^^ [Beelen et al.1 ([20061 ). 



Signatures of SN-origin dust hav e been claimed in the 
extincti on curves of a z ^ 6 .2 QSO (jMaiolino et al.l [20041 : 
see also [Gallerani et "aI1[2010[ ) and of two gamma- ray burst 
host galaxies, one at z ~ 6.3 ([Strat ta ct al.| [2007l but this 
re sult was undermine d bv [Zafar et a l. 2010) and one at z ~ 
5 ([Perlev et al.[[2009[) . 

The objective of this paper is to investigate if SNc and 
AGB stars are efficient enough to form dust at rcdshifts 
5 < z < 6.5 (1.15-0.85 Gyr after the Big Bang), or if grain 
growth in the ISM is required. We use a cosmological model 
with iJo = 70 km s^^ Mpc-\ and (Ha, f7,„) = (0.7, 0.3). 



2. Methodology 

In order to constrain the dust production efhciency in the 
early Universe, we selected z > 5 QSOs detected in both the 
millimeter continuum and CO lines allowing estimates to 
be made of their dust, gas and dynamical masses (Tab.[T]). 

We calculated dust masses (Mdust) from the 
1200 /xm data (re s t-fram e 160-200 /im) using Eq. 5 of 
[Michalowski et al.l ([2009[ ) assuming /? = 1.3. For three 
QSOs we adopted the derived dust temperatures (Tdust; 
Tab. [T|). For the rest we assumed the average of these 
estimates Tjust = 50 K. We assumed the mass absorption 
coefficient Ki2nn//.m = 0.67 cm^g~^, a conservatively high 
value (cf. [Alton et al.[[2004 ) resulting in systematically low 

In order to explore the impact of systematic un- 
certainties on Md u st 1 we also assumed / 3 = 2.0 (see 
Dunne et al.l [2000[ : |D unne fc Baled [200l[: [Vlahakis et all 
20051) . This gi ves Mdust smaller by a factor of ~ 3.75 



(see Fig. 3 of [Michalowski et al.l [2010a[) . Changing Tdust 
to a very hig h value of 80 K (compare with Fig. 2 of 



OSOs ([Haas et all 19981: Benford et al. 19991 Leech et al.l 


2OOII: Priddey & McMahonI [2001; Knudsen 


L et al.l 20031: 


Beelen et al.[ [20061: [Aravena et al.[ 20081 


Leioski et al.l 


20101). decreases the Mri,,=t bv a factor of ^ 


^ 2.3. Hence, 



we also assumed (Tdust, /3) = (80, 2.0). This results in strict 
lower limits on Mdust smaller by a factor of 3. 75 x2. 3 = 8. 6. 
However, this very conservative assumption is only chosen 



to illustrate an extreme limit on Mdust- It is not likely that 
the real values arc close to this limit as Tdust has been con- 
strained to be below 60 K for four out of nine QSOs in 
our sample with good wavelength coverage in the infrared 
(Tab.d]). 

Similar to [Wang et al.[ (j2010t) . we assume that the stellar 
masses (M* ) of the QSO host galaxies can be approximated 
as the difference between the dynamical (Mdyn; i-C- total) 
and the H2 gas masses (Mgas). The true values of M* are 
lower, unless QSOs harbour very little atomic gas (HI). 
Given the significant uncertainties in the conversion from 
CO line strength to Mgas, we also performed the calcula- 
tions with the upper limit setting M* equal to Afdyn- 

The inclination an gle of the gas disk , i, was adopted 
to be 65° for QSO 6 ([Walter et al.[[2004D and 40° for the 
others ([Wang et al.|[2010t) . The latter assumption is a major 
source of uncertainty in our analysis and is discussed below. 

We calculated the dust yields per AGB star and per 
SN (amount of dust formed in ejecta of one star) required 
to explain the inferred dust mas s es in t he z > 5 QSOs 
as described in [Michalowski et al.l (j2010bl ). The number of 
stars with masses between Mq and Mi in the stellar pop- 
ulation with a total mass of Mh. was calculated as N{Mq- 
Mi) = JI^^ M-'^dM/ JI^^I^I" M-°'MdM. We adopted 
an IMF with Mmin = 15, M^ax = 120 Mg, and a slope 
a = 2.35 (|SalDeteil[l955i or a = 1.5 for a top-heavy IMF). 
The average dust yield per star is Mdust/A^(Mo-Mi). 



3. Results and Discussion 

First, we consider a single dust producer i.e. assume that 
dust in the z > 5 QSOs was produced by either AGB stars 
or SNe. The required dust yields per AGB star and per 
SN are listed in Tab. [5| and shown in Fig. [T] as a function 
of redshift. Circles correspond to reasonable estimates of 
Tdust, /? and M*, whereas other values are shown to quan- 
tify the impact of the systematic uncertainties (error bars 
extend down to the reasonable lower limits, whereas arrows 
represent strict and unlikely lower limits). 

Except for QSOs 2 and 8 the required yields for AGB 
stars exceed the theoretically allowed maximum values 
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Theoretical yields 




Top-heavy IMF 



Theoretical yields: 
without dust destruction 
■ with dust destruction 

- - Cassiopeia A and Kepler 
// Other SN remnants 




5.0 5.2 5.4 



5.6 5.8 
Redshift 



6.0 6.2 6.4 



Fig. 1. Dust yields per AGB star (top) or per SN {bottom) required to explain dust in the z > 5 QSOs. For reasonable 
assumptions on the dust properties, AGB stars are not efficient enough and SNe would need to be unfeasibly efficient to 
form dust in these sources suggesting rapid grain grown in the ISM is likely to be responsible for the large dust masses. 
Circles: the best estimates of the required dust yields with error bars reflecting the uncertainty of /3 and A/*. Numbers 
indicate the QSOs as in Tab. [1] Arrows: strict and unlikely lower limits with very high Tdust and shown where data 
allow it (Tab. [T]). Gray symbols indicate that a top-heavy IMF was adopted. Dashed line and diagonal lines: the dust 
yields derived for Cassiopeia A, Kepler (~ IMq) and other SN remnants ('^ 10~^-10~^ Mq), respectively. Green area: 
theoretical dust yields for AGB stars (< 4 • 10~^ Mq). Light blue and blue areas: theoretical SN dust yields without 
(< 1.32 Mq) and with dust destruction implemented (< 0.1 Mq), respectively. 



(green area) by a factor of 2-15. The yields remain too 
high even for M^, = Mdyn and (5 = 2. They are consis- 
tent (though at the high end) with the theoretical predic- 
tions only under the unrealistic assumption of (Tdust,/?) = 
(80,2.0). Using the Tdust limits (Tab. [1]), we can robustly 
rule out a significant contribution of AGB stars to the dust 
formation in five out of nine QSOs (1, 3, 4, 5 and 6) and 
rule out their contribution in QSOs 7 and 9, unless their 
emission is dominated by hot (~ 80 K) dust. 

Therefore AGB stars are not efficient enough to form 
dust in the majority of the z > 5 QSOs. This contradicts 
the claim of IVahante et all (12009, ') that - 80% of dust in 
QSO 6 was created by AGB stars. The disagreement can 
be traced to the fact that they assumed M* ^ 10^^ M0, 
exceeding the Mdyn by a factor of ~ 15. 

For only two QSOs (2 and 8) are the required SN dust 
yields marginally within the theoretically predicted lim- 
its with dust destruction implemented (dark blue area on 
Fig. [U. For the remaining seven QSOs, one would need to 
assume unrealistically high Tdust md steep spectral slopes 



and in some cases an IMF more top-heavy than the Salpeter 
IMF. 

For these seven QSOs (including QSO 5 for which 
iMaiolino et"alll2004l claimed SN -origin dust) the required 
SN dust yields are within the theoretical limits without 
dust destruction (light blue area) and the values observed 
for SN remnants Cassiopeia A and Kepler (dashed line) . 

We checked that allowing AGB stars to form only a 
fraction of dust in the z > 5 QSOs and assigning the rest 
to SNe may have an impact on our conclusions for only 
four out of nine QSOs (3, 5, 7 and 9). This is illustrated 
in Fig. [21 where we show the required dust yields assum- 
ing different fractions of dust attributed to SNe. Solid lines 
represents the required yields for the QSOs. An increase 
in the fraction of SN dust corresponds to moving towards 
bottom-right corner (i.e. higher SN yields and lower AGB 
yields are required). If a curve corresponding to a QSO 
crosses the hatched region, corresponding to the allowed 
yields for both AGB stars and SNe, then these stellar ob- 
jects can account for dust in this QSO. Hence we conclude 
similarly as before, that combined AGB stars and SNe are 
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Table 2. Dust yields per star required to explain dust in 2; > 5 QSOs 



Td Dust Yields [Mq Per StarJ 



Dust Producer 


(K) 


P 


IMF 






1 


2 


3 


4 


5 


6 


7 


8 


9 


Sym 


AGB (2.5-8Af0) 


50 


1.3 


Sal. 


Mdyn 


-Mgas 


0.42 


0.02 


0.08 


0.64 


0.13 


0.45 


0.18 


0.03 


0.11 


• 


AGB (2.5-8M0) 


50 


2.0 


Sal. 


Mdyn 




0.08 


0.01 


0.02 


0.11 


0.03 


0.08 


0.03 


0.01 


0.03 


_L 


AGB (2.5-8Mq) 


80 


2.0 


Sal. 


Mdyn 




0.031 


0.003 


0.009 


0.048 


0.012 


0.039 


0.015 


0.003 


0.011 


t 


SN (8-4OM0) 


50 


1.3 


Sal. 


Mdyn 


-Mgas 


1.79 


0.11 


0.34 


2.76 


0.55 


1.93 


0.76 


0.12 


0.48 


• 


SN (8-4OM0) 


50 


1.3 


Top 


Mdyn 


-Mgas 


0.76 


0.05 


0.15 


1.18 


0.23 


0.82 


0.32 


0.05 


0.21 


• 


SN (8-4OM0) 


50 


2.0 


Sal. 


Mdyn 




0.35 


0.03 


0.08 


0.46 


0.12 


0.34 


0.15 


0.03 


0.11 


_L 


SN (8-4OM0) 


80 


2.0 


Sal. 


Mdyn 




0.13 


0.01 


0.04 


0.21 


0.05 


0.17 


0.06 


0.01 


0.05 


t 


SN (8-4OM0) 


80 


2.0 


Top 


Mdyn 




0.057 


0.005 


0.016 


0.088 


0.023 


0.072 


0.028 


0.005 


0.020 


t 



Notes. The IMF is either ISalpeteil l|l955f ) with a = 2.35 or top-heavy with a = 1.5. The M« column indicates either that stellar 
mass was assumed to be the difference between the dynamical and gas masses (Mdyn — Mgas) or that the strict upper limit to the 
stellar mass equal to the dynamical mass was adopted (see Sec. [2]). The numbered columns contain the required dust yields for all 
QSOs in the order given in Table [1] Only their numbers are given for brevity. The last column gives the symbol used on Fig. [T] 



efficient enough to form dust in QSOs 2 and 8 and are not 
efficient enougli for QSO 1, 4, and 6. The stellar dust pro- 
ducers may account for dust in QSOs 3, 5, 7 and 9, but only 
if very little dust is destroyed in SN shocks (light blue area 
on Fig. [2]). For these cases, SNe should be responsible for 
more than 50-75% of dust in these QSOs. Alternatively, 
the required dust yields for these four QSOs can be rec- 
onciled with theoretical expectations with dust destruction 
implemented (dark blue area on Fig.[2|) if wc assume a high 
value oi j3 — 2. 

We stress that our results are sensitive to the assumed 
gas disk inclinations. The required AGB and SN dust yields 
for individual QSOs decrease to theoretically allowed values 
(with dust destruction) for inclinations lower than 5-20°. 
It is however u nlikely that all our Q SOs exhibit such low 
inchnation (e.g. iPoUetta et al.ll2008l did not find any pre- 
ferred inclination for luminous QSOs). At least this is not 
the case for QSO 6 with a measured inclination of ~ 65°. 

Moreover, our derived required dust yields should be 
corrected towards lower values if i) the gas disk radius of 
QSOs is larger than 2.5 kpc assumed bv lWang et al.l (|2010l 
then the dynamical mass would be larger), or ii) the stellar 
component is more extended than the gas disk (then our 
upper limit on equal to Mdyn would apply only to the 
stellar component distributed within the extent of the gas 
disk). 

It is however unlikely that these conditions are ful- 
filled in our sample. Using the high-resolution CO line 
observations, the sizes of the gas disks have bee n con- 
strained for QSO 1 (< 3 kpc : iMaiolino et al.l[2007t ). QSO 
5 (2.2 X 5.0 kpc; 



Wang et al.ll2010D and OSO 6 (2.5 kpc 



IWalter et al.ll2004 ). Moreover, the star- forming gas of QSO 
6 ha s been found to be distributed within a radius of 0.75 
kpc (jWalter et al.ll2009f) . 

There is no estimate of the extent of the stellar com- 
ponent of the z > 5 QSOs, but at redshifts ^0 — 3 QSOs 
are typically hosted in 5, 3 kpc galaxies (jRidgwav et al.l 



I2OOII: IVeilleux et ani2009l). consist ent with a value of 2.5 
kpc assumed bv lWang' "et all (|2010[ ). 

Hence, we conclude that, unless the inclinations are 
biased low or the extent of stellar component are signifi- 
cantly larger than 2.5 kpc, both AGB stars and SNe would 
have to form unfeasibly large amounts of dust to account 
for dust present in the z > 5 QSOs. This may be taken 
as an indication of another (non-stellar) dust source in 
these objects, e.g., significant grain growth in the ISM 



(e.g^ lDraine fc SalDeteij|1979t iDwek fc ScalolllQSOl: iDraind 
120091 1 on the dust seeds produced by SNe or possibly AGB 
stars. Assuming that star formation in these QSOs began 
at z ~ 10, a timescale for in situ grain growth of a few x 10 
Myr jDrainelll990l . [20091 : iHirashital I2OOOI: IZhukovska et al.l 
I2OO8D is ^10% of the available time, suggesting ample time 
for grain growth in the ISM to explain the observed dust 
masses. 

Do stellar sources deliver enough additional heavy el- 
ements (not incorporated in dust) necessary for grain 
growth? The majority of heavy elements produced by 
AGB stars are already bound in dust grains (yields 
of carbon and other heavy e lemen ts are < 3.5 x 
10-2 Mq; iMorgan fc Edmunds! I2003D . On the other 
hand, a S N produces as m uch as < 1 M^ of heavy 
elements (iTodin i fc Fcrra ral 1200 it iNozawa et all 120031: 
iBianchi fc Schneider.,2007i :' lCherchneff fc Dwekll2009D . close 
to the required yields for the z > 5 QSOs (lower panel of 
Fig. [1]). Hence, even though SNe themselves do not pro- 
duce enough dust, they may deliver enough heavy elements 
to fuel the dust grain growth in the ISM. 



4. Conclusions 

We have derived the dust yields per AGB star and per SN 
required to explain observationally determined dust masses 
in 5 < 2; < 6.5 QSOs. We find that the yields for AGB stars 
typically exceed the theoretically allowed values making 
these objects inefficient to produce dust at high redshifts. 
SNe could in principle be responsible for dust in some of 
the QSOs, but with a requirement of high dust yields. This 
advocates for non-stellar dust source e.g. significant dust 
grain growth in the ISM of at least three out of nine QSOs. 
Wc argue that SNe deliver enough heavy elements to fuel 
the dust growth. 
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per. The Dark Cosmology Centre is funded by the Danish National 
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